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Введение. Множество практических задач опирается на задачу 
покрытия множеств: построение расписаний, расположение 
пунктов обслуживания, построение электронных схем. Это 
определяет актуальность поиска способов повышения 
эффективности решения данной задачи. 

Материалы и методы. Рассматриваются методы решения 
задачи о покрытии множества точным и приближенным 
алгоритмами. В качестве приближенного метода используется 
генетический алгоритм, в качестве точного — метод ветвей и 
границ. 

Результаты исследования. Генетический алгоритм во всех 
своих модификациях по временным характеристикам показал 
предсказуемость и стабильность во всех сериях 
экспериментов. Метод ветвей и границ был применен к задаче 
покрытия множеств и показал точные результаты. 
Обсуждение и заключения. Проведенные исследования 
показали, что для множеств небольших размеров 
целесообразно использовать метод ветвей и границ, который 
продемонстрировал быстрое время выполнения при 
гарантированно точном результате. Для множеств больших 
размеров рекомендуется использовать — генетический 
алгоритм, который гарантирует результат с незначительной 
погрешностью, причем изменение времени его работы 
стабильно и предсказуемо. 
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Введение. Наиболее важным классом задач в теории алгоритмов является класс М№Р-полных задач. На 
сегодняшний день нахождение оптимальных решений в этой области не перестает быть актуальным. В частности, 
большое количество прикладных задач можно свести к задаче о покрытии множества. Примерами могут являться 
многие задачи дискретной оптимизации: задачи стандартизации, упаковки и разбиения множества, задача о 
наибольшей клике, построение расписаний. Известна также и обратная сводимость задачи о покрытии к этим задачам. 





* Работа выполнена в рамках инициативной НИР. 
"" Е-пай: е-та]:х1еогх92(@тай.ти, с}-х@уапдех та, уа]ега33305@тай.та 
"Те тезеагсь 15 допе И т Фе Нате оЁ фе п14ерепаети В &0. 


Информатика, вычислительная техника и управление 


— 
9х) 
1 


//уезщК.допзва.га 


Вр 


138 


Вестник Донского государственного технического университета 2017, №3(90), 137-144 





На практике задачи о покрытии возникают при размещении пунктов обслуживания, в системах 
информационного поиска, при назначении экипажей на транспорте, проектировании интегральных схем и 
конвейерных линий ит. д. 


Постановка задачи. Дано множество И из и элементов и набор подмножеств (И, 5={$1,..., эй. Каждому 
подмножеству 5; сопоставлена некоторая неотрицательная стоимость с: 5—0*. 5'<5 является покрытием, если любой 
элемент из И принадлежит хотя бы одному элементу из 5" [1-3]. 

Задача о покрытии множествами заключается в нахождении набора подмножеств, покрывающего все 
множество И и имеющего минимально возможный вес (в случае взвешенной задачи) или минимально возможное 
число подмножеств (в случае невзвешенной задачи). 

Можно представить задачу в матричном виде [4, 5]. Пусть А=(а,) — произвольная матрица размера тхи с 
элементами а,Е {0,1} без нулевых строк и столбцов. Будем говорить, что в А строка { покрывается столбцом /, если 
а;=1. Подмножество столбцов называется покрытием, если в совокупности они покрывают все строки матрицы А. 
Пусть каждому столбцу поставлено в соответствие положительное число с, называемое весом столбца. Требуется 
найти покрытие минимального суммарного веса. Вводя переменные х, равные единице, если столбец ] входит в 
искомое покрытие, и равные нулю в противном случае, приходим к следующей формулировке задачи о покрытии: 

п 
хе „х, > шш, (1) 
при ограничениях 
я 
а (2) 
те 

Целью данного исследования является выявление целесообразности использования некоторых алгоритмов 
для решения задачи о покрытии множества на основе сравнения эффективности их работы по критериям 
быстродействия и точности. Для сравнения были выбраны: генетический алгоритм, метод ветвей и границ, алгоритм 
полного перебора. 


Генетический алгоритм. Генетические алгоритмы основаны на идеях естественного отбора и наследования 
и относятся к вероятностным эвристическим методам [6-9]. 

В настоящем исследовании используется модифицированная модель Голдберга: оптимизационной функцией 
является минимизация веса покрытия, а условием останова — неизменность лучшего решения в течение заданного 
числа поколений. Особенностью исследуемого алгоритма является то, что в кроссовере участвуют все особи текущего 
поколения. Кроме того, для каждого потомка применяется оператор мутации. 

Рассмотрим механизм кодирования особи. Каждый индивид К представлен хромосомой, являющейся п- 
мерным вектором хь, у которого /-й элемент ху, принимает значение 1, если подмножество 5, входит в покрытие, а в 
ином случае принимает значение 0. С таким представлением степень приспособленности }. индивида х, может быть 
рассчитана следующим образом: 

п 
л=Ус р (3) 
= 
где с‚— стоимость подмножества 5, 

Таким образом, оптимизационная функция выглядит как /— 1111. 

Для выбора родительских особей используется случайный отбор. Также применяется стратегия элитного 
отбора для улучшения точности результата без временных затрат: элитная особь формируется в первом поколении с 
помощью жадного алгоритма, описанного в статье [1]. 

В алгоритме используется тип мутации, основанный на изменении случайного гена на противоположное 
значение. Оператор скрещивания точечный. Выбираются пары хромосом из родительской популяции. Далее для 
каждой пары отобранных таким образом родителей разыгрывается позиция гена (локус) в хромосоме, определяющая 
так называемую точку скрещивания — [,. В результате скрещивания пары родительских хромосом получается 
следующая пара потомков: Р‚ — потомок, хромосома которого на позициях от | до [, состоит из генов первого 
родителя, а на позициях от [+1 до Г, — из генов второго родителя; Р› — потомок, хромосома которого на позициях от 
1 до [+ состоит из генов второго родителя, а на позициях от [+1 до Г — из генов первого родителя. Начальное 
поколение формируется из особей, соответствующих найденным случайным образом покрытий. 

Особенность работы генетического алгоритма состоит в том, что при скрещивании и мутации могут 
появиться особи, соответствующие покрытия которых не существуют, то есть являющиеся недопустимыми 
решениями. От этих решений необходимо избавляться. Алгоритм проверяет, существует ли покрытие, и если нет, то 
пытается, в случае скрещивания, выбрать другую вторую родительскую особь, а в случае мутации — выбрать другой 
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ген для его инвертирования. Если и это не «исправит» особь, родитель выберется заново случайным образом. 
Потомок заменяет случайно выбранную особь, если его приспособленность выше. 


Алгоритм полного перебора (брутфорс). Полный перебор — точный метод решения оптимизационных 
задач, относящийся к классу методов поиска решения исчерпыванием всевозможных вариантов. Сложность полного 
перебора зависит от количества всех возможных решений задачи. 

Полным перебором можно решить любую задачу из класса МР-полных задач. В зависимости от количества 
всех возможных решений задачи и времени вычисления целевой функции от каждого решения полный перебор может 
потребовать экспоненциального времени работы. 

Идея алгоритма полного перебора для решения задачи покрытия заключается в следующем: 

1) поиск всех возможных сочетаний подмножеств исходного множества и сравнение их целевых функций; 

2) определение покрытий среди найденных сочетаний; 

3) нахождение покрытия минимального веса. 

Из теории множеств известно, что число всех подмножеств множества из и элементов равно 2”. Другими 
словами, 


УС =”. (4) 
К=0 


Метод ветвей и границ. Данный метод является модификацией алгоритма полного перебора, что 
гарантирует точность результата его работы. Суть метода заключается в построении дерева полного перебора и 
отсечении бесперспективных ветвей решения по мере его обхода, что существенно уменьшает время его работы [10]. 
Условия задачи описываются формулами (1) и (2). Дерево перебора в этом случае строится так, что на каждом уровне 
К для каждого из узлов К-1 уровня добавляются в качестве дочерних узлов все возможные варианты ар-х,. Для 
поставленной задачи на каждом уровне добавляется количество узлов, равное количеству а,70. Алгоритм полного 
перебора предполагает полный обход такого дерева. Алгоритм после обработки каждого узла А-ого уровня приступает 
к обработке дочерних узлов прежде, чем переходит к следующему узлу К-ого уровня. Метод ветвей и границ 
сокращает время поиска оптимального решения за счет того, что при вхождении в каждый узел выполняется верхняя 
и/или нижняя оценка возможного решения, к которому приведет обход поддерева, корнем которого является текущий 
узел. В соответствии с полученной оценкой делается вывод — если лучшее из возможных решений хуже текущего, то 
данное поддерево (ветвь) отсекается и обход продолжается со следующего узла того же уровня, на котором было 
произведено отсечение. 

Таким образом, ключевым аспектом работы алгоритма по методу ветвей и границ является эффективная 
оценка поддерева, обязательным требованием к которой также является недопустимость потери точности. Поскольку 
алгоритм оценки поддерева опирается на текущее решение, инициализация начального решения оставляется на 
усмотрение разработчика. Возможны четыре варианта инициализации: 


1) заполнение решения максимальными значениями; 

2) выбор случайного покрытия; 

3) последовательная инициализация лучшими значениями (алгоритм Алексеева); 
4) нахождение решения приближенным алгоритмом. 


Поскольку многообразие и эффективность работы приближенных алгоритмов являются отдельным 
предметом исследования, для данного исследования их использование в качестве начального значения недопустимо. 
В связи с этим был выбран вариант инициализации по алгоритму Алексеева: последовательно, на каждом шаге в 
решение добавляется такая переменная х,, которая дает на данном шаге лучшее решение. Если таких переменных 
несколько, выбирается первая найденная. 

Для задачи покрытия множеств был разработан следующий оценочный метод. Пусть М — текущее решение 
задачи, то есть минимальное найденное количество подмножеств 5, покрывающих исходное множество И. В таком 
случае, минимальное улучшение равняется М№-1. Текущим состоянием М назовем число подмножеств 5, 
покрывающих сформированный на данном этапе набор х1хо ...х;. При М=М-1 выполняется проверка: если среди 
непокрытых переменных х; есть хотя бы одна, которая не покрывается подмножествами 5, формирующими состояние 
М, то минимальное улучшение невозможно, следовательно, данная ветвь считается бесперспективной и 
отбрасывается. В противном случае текущее состояние М становится текущим решением задачи. 


Входные данные и исходные параметры алгоритмов. В качестве входных данных для всех алгоритмов 
взята матрица А, представленная формулой (2). 

Для сравнения эффективности работы были выбраны генетический алгоритм и метод ветвей и границ. 

В качестве исходных параметров генетического алгоритма установлены: вероятность скрещивания — 100%, 
вероятность мутации — 100%. Для исследования выбрано 3 варианта с различным количеством особей — 100, 200 и 
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300. Для улучшения точности результата используется принцип элитного отбора. Для метода ветвей и границ 
начальное решение формируется алгоритмом Алексеева. 


Подтверждение точности метода ветвей и границ. С целью подтверждения точности алгоритма по методу 
ветвей и границ, а также доказательства его эффективности по временным показателям, были проведены серии 
экспериментов по сравнению его работы с работой алгоритма полного перебора на матрицах различного размера. 

В представленных ниже таблицах приведены усредненное значение результата и среднее время работы в 
миллисекундах, алгоритм полного перебора обозначен «АПП», метод ветвей и границ — «ВиГ». 










































































Таблица 1 
ТаЫе 1 
Результаты экспериментов при размере матрицы А — 10х10; количестве экспериментов — 500 
Кезий; орехрентепт5 жий А тайлх 512е — 10х10; питВег оГехреттепи; — 500 
АПП ВиГ 
Результат 2,288 2,288 
Время, мс 0,8343 0,0181 
Таблица 2 
ТаЫе 2 
Результаты экспериментов при размере матрицы А — 15х15; количестве экспериментов — 500 
Кезий; орехрейтет жий А тах я12е — 15х15; питВег оГехреттет5 — 500 
АПП ВиГ 
Результат 2,398 2,398 
Время, мс 56,6256 0,0617 
Таблица 3 
ТаЫе 3 
Результаты экспериментов при размере матрицы А — 20х20; количестве экспериментов — 500 
Кезий; орехрентепт5 ий А тах 512е — 20х20; питВег оГехреттепт5 — 500 
АПП ВиГ 
Результат 2,612 2,612 
Время, мс 3133,0340 0,1183 
Таблица 4 
ТаЫе 4 


Результаты экспериментов при размере матрицы А — 25х25; количестве экспериментов — 500 
Кезий; орехрентет5 ий А тах 512е — 25х25; питЬег оГехреттет; — 500 




















АПП ВиГ 
Результат 2,614 2,614 
Время, мс 155340,0994 0,2492 








Четыре серии экспериментов (суммарно 1600) не выявили ни одного расхождения результатов работы 
алгоритмов полного перебора и метода ветвей и границ. Таким образом, точность метода ветвей и границ доказана не 
только теоретически, но и экспериментально. 

По временным показателям метод ветвей и границ сработал быстрее во всех 4 сериях экспериментов: для 
матрицы размером 10х10 — в 46 раз, 15х15 — в 918 раз, 20х20 — в 26 480 раз и для 25х25 — в 623 451 раз. 

Данные результаты позволяют сделать однозначный вывод о преимуществе использования метода ветвей и 
границ перед алгоритмом полного перебора и в дальнейшем в исследовании в качестве точного алгоритма будет 
использоваться только метод ветвей и границ. 


Коновалов И. С. и др. Сравнение эффективности работы точных и приближенных алгоритмов 





Сравнение эффективности работы точных и приближенных алгоритмов. В данном исследовании для 
сравнения было выбрано 4 алгоритма: 3 модификации генетического алгоритма с использованием элиты, с 
количествами особей 100 (ГА100), 200 (ГА200) и 300 (ГАЗ00), а также алгоритм по методу ветвей и границ (ВиГ). 

В представленных ниже таблицах приведены усредненное значение результата (для вычисления среднего 


отклонения) и среднее время работы в миллисекундах. 


























































































































Таблица 5 
ТаЫе 5 
Результаты экспериментов при размере матрицы А — 20х20; количестве экспериментов — 500 
Кезий; орехрентет5 жий А тах 512е — 20х20; питЬег оГехреттет; — 500 
ГА1О0 ГА200 ГАЗ00 ВиГ 
Результат 2,688 2,684 2,684 2,684 
Время, мс 426,154 1033,16 1435,32 0,18 
Таблица 6 
ТаЫе 6 
Результаты экспериментов при размере матрицы А — 50х50; количестве экспериментов — 500 
Кезий; оГехрентет тий А тах 512е — 50х50; питЬег оГехреттепти5 — 500 
ГА1О0 ГА200 ГАЗ00 ВиГ 
Результат 3,112 3,024 3 3 
Время, мс 1406,05 3443,07 5266,25 8,21 
Таблица 7 
ТаЫе 7 
Результаты экспериментов при размере матрицы А — 75х75; количестве экспериментов — 500 
Кезий; оГехрейтет тий А тах я2е — 75х75; питЬег оГехреттепти5 — 500 
ГА1О0 ГА200 ГАЗ00 ВиГ 
Результат 3,846 3,684 3,572 3,142 
Время, мс 1724,13 3553,52 5437,28 320,64 
Таблица 8 
ТаЫе 8 
Результаты экспериментов при размере матрицы А — 85х85; количестве экспериментов — 200 
Кезий; оГехрейтепт5 лий А тах 512е — 85х85; питЬег оГехреттет; — 200 
ГА1о0 ГА200 ГАЗ00 ВиГ 
Результат 4,01 3,925 3,85 3,445 
Время, мс 2030,56 4159,29 6354,38 8851,48 
Таблица 9 
ТаЫе 9 


Результаты экспериментов при размере матрицы А — 85х85; количестве экспериментов — 400 


Кезий; оГехрентепт5 тий А тах 512е — 85х85; питЬег оГехреттет; — 400 


























ГА100 ГА200 ГАЗ00 ВиГ 
Результат 3,904 3,8275 3,8225 3,395 
Время, мс 2202,56 4312,41 6394,33 2892,97 
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Таблица 10 
ТаЫе 10 
Результаты экспериментов при размере матрицы А — 90х90; количестве экспериментов — 300 
Везий; оГехрентепт5 тий А тах 12е — 90х90; питЬег оГехреттепи; — 300 
ГА1О0 ГА200 ГАЗ00 ВиГ 
Результат 3,96 3,93 3,91 3,56 
Время, мс 2080,22 4735,22 7088,39 1926,04 
Таблица 11 
ТаМе 11 


Результаты экспериментов при размере матрицы А — 100х100; количестве экспериментов — 100 


Кезий$ оГехрентет5 жий А тайлх 512е — 100х100; питБег оГехрептепту — 100 


























ГА100 ГА200 ГАЗ00 ВиГ 
Результат 3,98 3,99 3,95 3,79 
Время, мс 4448,63 9450,40 14296,89 3655,89 








Таблица 12 
ТаЫе 12 


Результаты экспериментов при размере матрицы А — 100х100; количестве экспериментов — 200 
Кезий$ орехрентет5 жий А тайлх 512е — 100х100; питБег оГехреттептя — 200 











ГА100 ГА200 ГАЗ00 ВиГ 
Результат 4,015 3,995 3,924 3,85 
Время, мс 4374,55 9291,33 15378,67 38932,69 























Анализ полученных результатов. В сериях экспериментов 1, 2 и 3 точный алгоритм по методу ветвей и 
границ показал убедительное преимущество в скорости работы, превзойдя самую быструю модификацию 
генетического алгоритма (ГА100) в 2367,5 раз в первом случае, в 171,26 во втором ив 52,82 в третьем. 

Однако, начиная с четвертой серии экспериментов, время работы алгоритма по методу ветвей и границ 
сильно варьируется и показывает весь спектр результатов от лучшего к худшему. Объясняется это тем, что ход работы 
алгоритма сильно зависит от входной матрицы. При большой плотности единиц разветвленность дерева больше, а 
также больше ветвей, в которых оценка возможна лишь ближе к листьям. То есть, несмотря на то, что время работы 
метода ветвей и границ будет всегда меньше, чем время работы алгоритма полного перебора Т, в худшем случае 


Т 
время данного метода может стремиться к = › где и — некоторая константа. 


Генетический алгоритм во всех модификациях по временным характеристикам показал предсказуемость и 
стабильность во всех сериях экспериментов. Среднее отклонение от точного результата для ГА1О0 составило 9,77%, 
для ГА200 — 8,13%, для ГАЗО0 — 6,68%. Вышеизложенные результаты позволяют сделать вывод о том, что с ростом 
количества особей улучшается точность генетического алгоритма. 


Выводы. Проведенные исследования показали, что для множеств небольших размеров (до 75 элементов 
включительно) целесообразно использовать метод ветвей и границ, который продемонстрировал высокое 
быстродействие при гарантии точного результата. Для множеств больших размеров рекомендуется использовать 
генетический алгоритм, который гарантирует получение результата с незначительной погрешностью, но за 
определенный фиксированный промежуток времени. Метод же ветвей и границ, по результатам проведенных 
экспериментов, с увеличением размера входных данных может вести себя непредсказуемо. 
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